All-sky search for gravitational-wave bursts in the first joint LIGO-GEO-Virgo run 
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We present results from an all- sky search for unmodeled gravitational-wave bursts in the data 
collected bv the LIGO, GEO 600 and Virgo detectors between November 2006 and October 2007. 
The search is performed by three different analysis algorithms over the frequency band 50 - 6000 Hz. 
Data are analyzed for times with at least two of the four LIGO-Virgo detectors in coincident oper- 
ation, with a total live time of 266 days. No events produced by the search algorithms survive the 
selection cuts. We set a frequentist upper limit on the rate of gravitational-wave bursts impinging 
on our network of detectors. When combined with the previous LIGO search of the data collected 
between November 2005 and November 2006, the upper limit on the rate of detectable gravitational- 
wave bursts in the 64-2048 Hz band is 2.0 events per year at 90% confidence. We also present event 
rate versus strength exclusion plots for several types of plausible burst waveforms. The sensitivity 
of the combined search is expressed in terms of the root-sum-squared strain amplitude for a variety 
of simulated waveforms and lies in the range 6 X 10~ 22 Hz -1 / 2 to 2 x 10 -20 Hz“ l ^ 2 . This is the first 
untriggered burst search to use data from the LIGO and Virgo detectors together, and the most 
sensitive untriggered burst search performed so far. 

PACS numbers: 04.80.Nn, 07.05.Kf, 95.30.Sf, 95.85.Sz 


I. INTRODUCTION 

The LIGO Scientific Collaboration (LSC) and the 
Virgo Collaboration operate a network of interferomet- 
ric gravitational-wave (GW) detectors with the goal of 
detecting gravitational waves from astrophysical sources. 
Some of these sources may produce transient “bursts” 
of GW radiation with relatively short duration (<1 s). 
Plausible burst sources [1] include merging compact bi- 
nary systems consisting of black holes and/or neutron 
stars [2, 3], core-collapse supernovae [4], neutron star col- 
lapse [5], starquakes associated with magnetar flares [6] 
or pulsar glitches [7], cosmic string cusps [8], and other 
violent events in the Universe. 

During the most recent datar-taking run five GW de- 
tectors were operational. The three LIGO detectors [9] 
started their Science Run 5 (S5) in November 2005, and 
the GEO 600 detector [10] joined the S5 run in January 
2006. The Virgo detector [11] began its Virgo Science 


Run 1 (VSR1) in May 2007. All five instruments took 
data together until the beginning of October 2007. 

An all-sky search for GW burst signals has already 
been conducted on the first calendar year of the LIGO 
S5 data (referred to as c S5yl”) in a wide frequency band 
of 64 - 6000 Hz [12, 13]. In this paper, we report on 
a search for GW burst signals in the frequency band 
50 — 6000 Hz for the rest of the S5/VSR1 run, referred to 
as “S5y2/VSR1”. It includes data collected by the LIGO 
and Virgo detectors, which had comparable sensitivities, 
and uses three different search algorithms. In compar- 
ison with the S5yl analysis, the network of LIGO and 
Virgo detectors, spread over three sites, provides better 
sky coverage as well as improved capabilities to reject 
spurious signals. S5y2/VSR1 is also the first long-term 
observation with the world-wide network of interferomet- 
ric detectors. This is a major step forward with respect 
to previous observations led by the network of resonant 
detectors [14, 15], since, as we will show in this paper, the 
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performance is improved by more than one order of mag- 
nitude both in the analyzed frequency bandwidth and 
the level of instrumental noise. 

This paper is organized as follows. In Section II we 
describe the LSC and Virgo instruments. In Section III 
we give a brief overview of the search procedure. In Sec- 
tion IV we present the search algorithms. Simulations 
are described in Section V, and the error analysis in Sec- 
tion VL The results of the search are presented in Sec- 
tion VII, and astrophysical implications are discussed in 
Section VIII. The appendices provide additional details 
on data characterization and the analysis pipelines. 


II. DETECTORS 

A. LIGO 

LIGO consists of three detectors at two observato- 
ries in the United States. Each detector is a large 
Michelson-type interferometer with additional mirrors 
forming Fabry-Perot cavities in the arms and a power- 
recycling mirror in the input beam path. Interferometric 
sensing and feedback is used to “lock” the mirror po- 
sitions and orientations to keep all of the optical cavi- 
ties on resonance. A gravitational wave is sensed as a 
quadrupolar strain, measured interferometricaily as an 
effective difference between the lengths of the two arms. 
The LIGO Hanford Observatory, in Washington, houses 
independent detectors with the arm lengths of 4 km and 
2 km, called HI and H2 respectively. The LIGO Liv- 
ingston Observatory, in Louisiana, has a single detector 
with 4-km arms, called LL The detector instrumenta- 
tion and operation are described in detail elsewhere [9], 
and the improvements leading up to the S5 run which are 
most relevant for GW burst searches have been described 
in the first-year search [12]. 

The best achieved sensitivities of the LIGO detectors 
during the second year of S5, as a function of signal fre- 
quency, are shown in Fig. 1. The detectors are most sen- 
sitive over a band extending from about 40 Hz to a few 
kHz. Seismic noise dominates at lower frequencies since 
the effectiveness of the seismic isolation system is a very 
strong function of frequency. Above ~200 Hz, laser shot 
noise corrected for the Fabry-Perot cavity response yields 
an effective strain noise that rises linearly with frequency. 
The sensitivity at intermediate frequencies is determined 
mainly by thermal noise, with contributions from other 
sources. The peaks at ~350Hz and harmonics are the 
thermally-excited vibrational modes of the wires from 
which the large mirrors are suspended. Smaller peaks 
are due to other mechanical resonances, power line har- 
monics, and calibration signals. 

Commissioning periods during the second year of S5 
led to incremental improvements in the detector sensi- 
tivities. The most significant of these were in January 
2007, when the seismic isolation systems at both sites 
were improved to reduce the coupling of microseismic 


noise to the mirror suspensions, thereby mitigating noise 
from the nonlinear Barkhausen effect [16] in the magnets 
used to control the mirror positions; and in August 2007, 
when the LI frequency stabilization servo was re-tuned. 
Overall, the average sensitivities of the HI and LI detec- 
tors during the second year were about 20% better than 
the first-year averages, while the H2 detector (less sensi- 
tive to begin with by a factor of ~2) had about the same 
average sensitivity in both years. The operational duty 
cycles for all three detectors also improved as the run 
progressed, from (72.8%, 76.7%, 61.0%) averaged over 
the first year to (84.0%, 80.6%, 73.6%) averaged over the 
second year for HI, H2, and LI, respectively. 


B. GEO 600 

The GEO 600 detector, located near Hannover, Ger- 
many, also operated during S5, though with a lower sensi- 
tivity than the LIGO and Virgo detectors. The GEO 600 
data are not used in the initial search stage of the current 
study as the modest gains in the sensitivity to GW signals 
would not offset the increased complexity of the analy- 
sis. The GEO 600 data are held in reserve, and used to 
follow up any detection candidates from the LIGO- Virgo 
analysis. 

GEO 600 began its participation in S5 on January 21 
2006, acquiring data during nights and weekends. Com- 
missioning work was performed during the daytime, fo- 
cussing on gaining a better understanding of the detector 
and improving data quality. GEO switched to full-time 
data taking from May 1 to October 6, 2006, then re- 
turned to night-and-weekend mode through the end of 
the S5 run. Overall GEO 600 collected about 415 days 
of science data during S5, for a duty cycle of 59.7% over 
the full S5 run. 


C. Virgo 

The Virgo detector [11], also called VI, is an inter- 
ferometer with 3 km arms located near Pisa in Italy. 
One of the main instrumental differences with respect 
to LIGO is the seismic isolation system based on super- 
attenuators [17], chains of passive attenuators capable of 
filtering seismic disturbances in 6 degrees of freedom with 
sub- Hertz corner frequencies. For VSR1, the Virgo duty 
cycle was 81% and the longest continuous period with the 
mirror positions interferometricaily controlled was more 
than 94 hours. Another benefit from super-attenuators is 
a significant reduction of the detector noise at very low 
frequency (< 40 Hz) where Virgo surpasses the LIGO 
sensitivity. 

Above 300 Hz, the spectral sensitivity achieved by 
Virgo during VSR1 is comparable to that of LIGO (see 
Figure 1). Above 500 Hz the Virgo sensitivity is domi- 
nated by shot noise. Below 500 Hz there is excess noise 
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FIG. i: Best noise amplitude spectral densities of the five 
LSC/Virgo detectors during S5/VSR1. 


due to environmental and instrumental noise sources, and 
below 300 Hz these produce burst-like transients. 

Due to the different orientation of its arms, the antenna 
pattern (angular sensitivity) of Virgo is complementary 
to that of the LI GO detectors, with highest response in 
directions of low LIGO sensitivity. Virgo therefore signif- 
icantly increases the sky coverage of the network. In ad- 
dition, simultaneous observations with the three LIGO- 
Virgo sites improve rejection of spurious signals and al- 
low reconstruction of the sky position and waveforms of 
detected GW sources. 


III. SEARCH OVERVIEW 

The analysis described in this paper uses data from 
the LIGO detectors collected from 14 November 2006 
through 1 October 2007 (S5y2), and Virgo data from 
VSR1, which started on 18 May 18 2007 and ended at 
the same time as S5 [40]. The procedure used for this 
S5y2/VSR1 search is the same as that used for S5yl [12]. 
In this section we briefly review the main stages of the 
analysis. 

A. Data quality flags 

The detectors are occasionally affected by instrumental 
or data acquisition artifacts as well as by periods of de- 
graded sensitivity or an excessive rate of transient noise 
due to environmental conditions such as bad weather. 
Low-quality data segments axe tagged with Data Quality 
Flags (DQFs). These DQFs are divided into three cate- 
gories depending on their seriousness. Category 1 DQFs 
are used to define the data segments processed by the 
analysis algorithms. Category 2 DQFs are unconditional 


data cuts applied to any events generated by the algo- 
rithms. Category 3 DQFs define the clean data set used 
to calculate upper limits on the GW rates. 

We define DQFs for S5y2/VSR1 following the ap- 
proach used for S5yl [12]. More details are gbnn in Ap- 
pendix A. After category 2 DQFs have been applied, the 
total available time during this period is 261.6 days for 
HI, 253.4 days for H2, 233.7 days for LI and 106.2 days 
for VI [41]. 


B. Candidate Event Generation 

As discussed in Section IV, three independent search 
algorithms are used to identify possible GW bursts: Ex- 
ponential Gaussian Correlator (EGC), fl-pipeline (H), 
and coherent WaveBurst (cWB). We analyze data from 
time intervals when at least two detectors were operat- 
ing in coincidence. Altogether, eight networks, or sets of 
detectors, operating during mutually exclusive time pe- 
riods are analyzed by at least one algorithm. Table I 
shows the time available for analysis (“live time”) for the 
different network configurations after application of cat- 
egory 1 and 2 DQFs. The actual times searched by each 
algorithm for each network (“observation times”) reflect 
details of the algorithms, such as the smallest analyz- 
able data block, as well as choices about which networks 
are most suitable for each algorithm. The three- and 
two-detector network configurations not shown in Table I 
have negligible live time and are not considered in this 
search. 


network 

live time 

cWB 

Cl 

EGC 

H1H2L1V1 

68.9 

68.2 

68.7 

66.6 

H1H2L1 

124.6 

123.2 

123.4 

16.5 

H1H2V1 

15.8 

15.7 

15.1 

15.3 

H1L1V1 

4.5 

4.2 

- 

4.4 

H1H2 

35.4 

35.2 

34.8 

- 

H1L1 

7.2 

5.9 

- 

- 

L1V1 

6.4 

- 

6.3 

- 

H2L1 

3.8 

3.5 

- 

- 


TABLE I: Exclusive live time in days for each detector net- 
work configuration after category 2 DQFs (second column) 
and the observation time analyzed by each of the search al- 
gorithms (last three columns). The cWB algorithm did not 
process the LI VI network because the coherent likelihood reg- 
ulator used in this analysis was suboptimal for two detectors 
with verv different orientations. Omega used a coherent com- 
bination of Hi and H2 as an effective detector and thus an- 
alyzed networks either with both or with neither. EGC ana- 
lyzed only data with three or more interferometers during the 
part of the run when Virgo was operational. 

LIGO and GEO 600 data are sampled at 16384 Hz, 
yielding a maximum bandwidth of 8192 Hz, while Virgo 
data are sampled at 20000 Hz. Because of the large cali- 
bration uncertainties at high frequency, only data below 
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6000 Hz are used in the search. Also, because of high seis- 
mic noise, the frequency band below 50 Hz is excluded 
from the analysis. Furthermore, the EGC search was 
limited to the 300-5000 Hz band over which Virgo’s sen- 
sitivity was comparable to LIGO’s. In Section VI we 
describe the influence of the calibration uncertainties on 
the results of the search. 


C. Vetoes 

After gravitational-wave candidate events are identi- 
fied by the search algorithms, they are subject to ad- 
ditional ‘Veto” conditions to exclude events occurring 
within certain time intervals. These vetoes are based 
on statistical correlations between transients in the GW 
channel (data stream) and the environmental and inter- 
ferometric auxiliary channels. 

We define vetoes for S5y2/VSR1 following the ap- 
proach used for S5yl [12]. More details are given in Ap- 
pendix B. 


D. Background Estimation and Tuning 

To estimate the significance of candidate GW events, 
and to optimize event selection cuts, we need to mea- 
sure the distribution of events due to background noise. 
With a multi-detector analysis one can create a sample 
of background noise events and study its statistical prop- 
erties. These samples are created by time-shifting data 
of one or more detectors with respect to the others by 
“un-physical” time delays ( i.e . much larger than the 
maximum time-of-flight of a GW signal between the de- 
tectors). Shifts are typically in the range from ~1 s to 
a few minutes. Any triggers that are coincident in the 
time-shifted data cannot be due to a true gravitational- 
wave signal; these coincidences therefore sample the noise 
background. Background estimation is done separately 
for each algorithm and network combination, using hun- 
dreds to thousands of shifts. To take into account pos- 
sible correlated noise transients in the HI and H2 detec- 
tors, which share a common environment and vacuum 
system, no time-shifts are introduced between these de- 
tectors for any network combination including another 
detector. 

The shifted and unshifted data are analyzed identi- 
cally. A portion of the background events are used to- 
gether with simulations (see below) to tune the search 
thresholds and selection cuts; the remainder is used to 
estimate the significance of any candidate events in the 
unshifted data after the final application of the selec- 
tion thresholds. All tuning is done purely on the time 
shifted data and simulations prior to examining the un- 
shifted data-set. This ‘blind” tuning avoid any biases in 
our candidate selection. The final event thresholds are 
determined by optimizing the detection efficiency of the 
algorithms at a fixed false alarm rate. 


E. Hardware and software injections 

At pseudo-random times during the run, simulated 
burst signals were injected (added) into the interferom- 
eters by sending pro-calculated waveforms to the mir- 
ror position control system. These “hardware injec- 
tions” provided an end-to-end verification of the detec- 
tor instrumentation, the data acquisition system and the 
data analysis software. The injection times were clearly 
marked in the data with a DQF. Most of hardware injec- 
tions were incoherent, i.e., performed into a single detec- 
tor with no coincident injection into the other detectors. 
Some injections were performed coherently by taking into 
account a simulated source location in the sky and the 
angle-dependent sensitivity of the detectors to the two 
wave polarization states. 

In addition to the flagged injections, a “blind injec- 
tion challenge” was undertaken in which a small num- 
ber (possibly zero) of coherent hardware injections were 
performed without being marked by a DQF. Informa- 
tion about these blind injections (including whether the 
number was nonzero) was hidden from the data analy- 
sis teams during the search, and revealed only afterward. 
This challenge was intended to test our data analysis pro- 
cedures and decision processes for evaluating any candi- 
date events that might be found by the search algorithms. 

To determine the sensitivity of our search to gravita- 
tional waves, and to guide the tuning of selection cuts, 
we repeatedly re-analyze the data with simulated signals 
injected in software. The same injections are analyzed 
by all three analysis pipelines. See Section V for more 
details. 


IV. SEARCH ALGORITHMS 

Anticipated sources of gravitational wave bursts are 
usually not understood well enough to generate wave- 
forms accurate and precise enough for matched filtering 
of generic signals. While some sources of GW bursts are 
being modeled with increasing success, the results tend to 
be highly dependent on physical parameters which may 
span a large parameter space. Indeed, some burst signals, 
such the white-noise burst from turbulent convection in 
a core-collapse supernova, are stochastic in nature and so 
are inherently not templatable. Therefore usually more 
robust excess-power algorithms [18-21] are employed in 
burst searches. By measuring power in the data as a 
function of time and frequency, one may identify regions 
where the power is not consistent with the anticipated 
fluctuations of detector noise. To distinguish environ- 
mental and instrumental transients from true GW sig- 
nals, a multi-detector analysis approach is normally used, 
in which the event must be seen in more than one detec- 
tor to be considered a candidate GW. 

The simplest multi-detector analysis strategy is to re- 
quire that the events identified in the individual detec- 
tors are coincident in time. The time coincidence win- 
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dow which should be chosen to take into account the 
possible time delays of a GW signal arriving at different 
sites, calibration and algorithmic timing biases, and pos- 
sible signal model dependencies. Time coincidence can 
be augmented by requiring also an overlap in frequency. 
One such time-frequency coincidence method used in this 
search is the EGC algorithm [22] (see also Appendix C). 
It estimates the signal-to-noise ratio (SNR) p* in each 
detector k and uses the combined SNR p C o m b = y/Yhk p\ 
to rank candidate events. 

A modification of the time-frequency coincidence ap- 
proach is used in the ft search algorithm [23] (also see 
Appendix D). In ft, the identification of the H1H2 net- 
work events is improved by coherently combining the HI 
and H2 data to form a single pseudo-detector data stream 
H_|_. This algorithm takes an advantage of the fact that 
the cc-located and co-aligned HI and H2 detectors have 
identical responses to a GW signal. The performance of 
the ft algorithm is further enhanced by requiring that no 
significant power is left in the H1-H2 null stream, H_, 
where GW signals cancel. This veto condition helps to 
reduce the false alarm rate due to random coincidences of 
noise transients, which typically leave significant power 
in the null stream. Network events identified by ft are 
characterized by the strength Z = p 2 / 2 of the individ- 
ual detector events, and by the correlated H1H2 energy 

r/corr 

z n+ - 

A different network analysis approach is used in 
the cWB search algorithm [24] (see also [12] and Ap- 
pendix E). The cWB algorithm performs a least-squares 
fit of a common GW signal to the data from the different 
detectors using the constrained likelihood method [25]. 
The results of the fit are estimates of the h+ and h x wave- 
forms. the most probable source location in the sky, and 
various likelihood statistics used in the cWB selection 
cuts. One of these is the maximum likelihood ratio L m , 
which is an estimator of the total SNR detected in the 
network. A part of the statistic depending on pair- 
wise combinations of the detectors is used to construct 
the network correlated amplitude 77, which measures the 
degree of correlation between the detectors. Random co- 
incidences of noise transients typically give low values 
of 17, making this statistic useful for background rejec- 
tion. The contribution of each detector to the total SNR 
is weighted depending on the variance of the noise and 
angular sensitivity of the detectors. The algorithm au- 
tomatically marginalizes a detector with either elevated 
noise or unfavorable antenna patterns, so that it does not 
limit the sensitivity of the network. 


V. SIMULATED SIGNALS AND EFFICIENCIES 

The detection efficiencies of the search algorithms de- 
pend on the network configuration, the selection cuts 
used in the analysis, and the GW morphologies which 
may span a wide range of signal durations, frequencies 
and amplitudes. To evaluate the sensitivity of the search 


and verify that the search algorithms do not have a strong 
model dependency, we use several sets of ad-hoc wave- 
forms. These include 

Sine-Gaussian waveforms: 


h+(t) - /iosin(2ir/oi) exp[— (2nfot) 2 /2Q 2 ], (5.1) 
MO = 0- (5.2) 

We use a discrete set of central frequencies / o from 
70 Hz to 6000 Hz and quality factors Q of 3, 9, and 
100; see Table II and Fig. 2 (top). The amplitude 
factor h 0 is varied to simulate GWs with different 
strain amplitudes. For definition of the polarize 
tions, see Eq. (5.8) and text below it. 

Gaussian waveforms: 

M*) = fcoexpM 2 /r 2 ), (5.3) 

M0 = 0, (5.4) 

where the duration parameter r is chosen to be one 
of (0.1, 1.0, 2.5, 4.0) ms; see Fig. 2 (middle). 

Harmonic ringdown signals: 

h+[t) = h 0> + cos(27r/ot) exp[— t/r], t > l/(4/ 0 jj5.5) 
M0 = Vx sin(27r/oOexp[-£/r], t> 0. (5.6) 

We use several central frequencies /o from 1590 Hz 
to 3067 Hz, one long decay time, r = 200 ms, 
and two short decay times, 1 ms and 0.65 ms; see 
Table III and Fig. 2 (bottom). Two polarization 
states are used: circular (ho,+ = ho,x)> and linear 
(/io,+ = 0). The quarter-cycle delay in h+ is to 
avoid starting the waveform with a large jump. 

Band-limited white noise signals: 

These are bursts of Gaussian noise which are white 
over a frequency band [/iow>/iow + A/] and which 
have a Gaussian time profile with standard devia/- 
tion decay time r; see Table IV. These signals are 
unpolarized in the sense that the two polarizations 
h+ and h x have equal RMS amplitudes and are 
uncorrelated with each other. 


The strengths of the ad hoc waveform injections are char- 
acterized by the root-square-sum amplitude h ra8 , 




dt(IM0l 2 + |M0l 2 )- (5-7) 


The parameters of these waveforms are selected to 
coarsely cover the frequency range of the search from 
~50 Hz to ~6 kHz, and duration of signals up to a few 
hundreds of milliseconds. The Gaussian, sine-Gaussian 
and ringdown waveforms' explore the space of GW signals 
with small time-frequency volume, while the white noise 
bursts explore the space of GW signals with relatively 
large time-frequency volume. Although the simulated 
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waveforms are not physical, they may be similar to some 
waveforms produced by astrophysical sources. For ex- 
ample, the sine-Gaussian waveforms with few cycles are 
qualitatively similar to signals produced by the mergers 
of two black holes [2]. The long- timescale ringdowns are 
similar to signals predicted for excitation of neutron-star 
fundamental modes [26] . Some stellar collapse and core- 
collapse supernova models predict signals that resemble 
short ringdown waveforms (in the case of a rapidly rotat- 
ing progenitor star) or band-limited white-noise wave- 
forms with random polarizations. In the context of the 
recently proposed acoustic mechanism for core-collapse 
supernova explosions, quasi-periodic signals of >500 ms 
duration have been proposed [4]. 

To test the range for detection of gravitational waves 
from neutron star collapse, two waveforms were taken 
from simulations by Baiotti et al. [5], who modeled neu- 
tron star gravitational collapse to a black hole and the 
subsequent ringdown of the black hole using collapsing 
polj tropes deformed by rotation. The models whose 
waveform we chose were Dl, a nearly spherical 1.67 A/ 0 
neutron star, and D4, a 1.86 M 0 neutron star that is 
maximally deformed at the time of its collapse into a 
black hole. These two specific waveforms represent the 
extremes of the parameter space in mass and spin con- 
sidered in [5]. They are linearly polarized (h x = 0), 
with the waveform amplitude varying with the inclina- 
tion angle i (between the wave propagation vector and 
symmetry axis of the source) as sin 2 6. 

The simulated detector responses hdet are constructed 
as 


/idet = F+(0, 4>, ij>)h+ + Fx {6, (5.8) 

Here F+ and F x are the detector antenna patterns, which 
depend on the direction to the source (0,0) and the 
polarization angle 0. (The latter is defined as in Ap- 
pendix B of [18].) These parameters are chosen ran- 
domly for each injection. The sky direction is isotrop- 
ically distributed, and the random polarization angle is 
uniformly distributed on [0,7r). The injections are dis- 
tributed uniformly in time across the S5y2/VSR1 run, 
with an average separation of 100 s. Note that for the 
ad-hoc waveforms no i is used. 

The detection efficiency after application of all selec- 
tion cuts was determined for each waveform type. All 
waveforms were evaluated using cWB, while subsets were 
evaluated using H and EGC, due mainly to the limited 
frequency bands covered by those algorithms as they were 
used in this search (48-2048 Hz and 300-5000 Hz, respec- 
tively;. Figure 2 shows the combined efficiency curves for 
selected sine-Gaussian, Gaussian and ringdown simulated 
signals as a function of the h TSS amplitude. Figure 3 shows 
the detection efficiency for the astrophysical signals Dl 
and D4 as a function of the distance to the source. 

Each efficiency curve is fitted with an empirical func- 
tion and the injection amplitude for which that function 
equals 50% is determined. This quantity, is a con- 
venient characterization of the sensitivity of the search to 
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FIG. 2: Efficiency for selected waveforms as a function of 
signal amplitude h raB for the logical OR of the H1H2L1V1, 
H1H2L1, and H1H2 networks. Top: sine-Gaussians with Q = 
9 and central frequency spanning between 70 and 5000 Hz. 
Middle: Gaussians with r between 0.1 and 4.0 ms. Bottom: 
linearly (L) and circularly (C) polarized ringdowns with t 
200 ms and frequencies between 1590 and 2590 Hz. 
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FIG. 3: Efficiency of the H1H2L1V1 network as a function 
of distance for the D1 and D4 waveforms of Baiotti et al [5] 
predicted by polytropic general-relativistic models of neutron 
star collapse. These efficiencies assume random sky location, 
polarization and inclination angle. 


that waveform morphology. Tables II, III, and IV sum- 
marize the sensitivity of the search to the sine-Gaussian, 
ringdown, and band-limited white noise burst signals. 
Where possible, we also calculate the sensitivity of the 
logical OR of the cWB and H algorithms (since those 
two are used for the upper limit calculation as described 
in Sec. VII) , and for the appropriately weighted combi- 
nation of all networks (some of which are less sensitive) 
contributing to the total observation time. In general, the 
efficiency of the combination of the search algorithms is 
slightly more sensitive than the individual algorithms. 

VI. UNCERTAINTIES 

The amplitude sensitivities presented in this paper, i.e. 
the hjss values at 50% and 90% efficiency, have been ad- 
justed upward to conservatively reflect statistical and sys- 
tematic uncertainties. The statistical uncertainty arises 
from the limited number of simulated signals used in the 
efficiency curve fit, and is typically a few percent. The 
dominant source of systematic uncertainty comes from 
the amplitude calibration: the single detector amplitude 
calibration uncertainties is typically of order 10%. Ne- 
glegibie effects are due to phase and timing uncertainties. 

The amplitude calibration of the interferometers is less 
accurate at high frequencies than at low frequencies, and 
therefore two different approaches to handling calibra^ 
tion uncertainties are used in the S5y2/VSR1 search. In 
the frequency band below 2 kHz, we use the procedure 
established for S5yl [13]. We combine the amplitude un- 
certainties from each interferometer into a single uncer- 
tainty by calculating a combined root-sum-square ampli- 
tude SNR and propagating the individual uncertainties 
assuming each error is independent: as a conservative re- 
sult, the detection efficiencies are rigidly shifted towards 
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[Hz] 

Q 
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9 

36.5 

• 

57.8 

36.5 

38.3 

146 

2000 

3 

42.6 

• 

- 

42.6 

47.1 

191 

2000 

9 

40.6 

- 

58.7 

40.6 

44.0 

177 

2000 

100 

34.9 

- 

- 

34.9 

38.4 

153 

2226 

9 

46.0 

- 

68.6 

46.0 

51.1 

187 

2477 

3 

61.9 

- 

- 

61.9 

65.6 

262 

2477 

9 

53.5 

• 

76.7 

53.5 

56.1 

206 

2477 

100 

44.5 

- 

- 

44.5 

48.9 

201 

2756 

9 

60.2 

- 

82.2 

60.2 

64.4 

248 

3067 

3 

86.9 

- 

- 

86.9 

87.0 

343 

3067 

9 

69.0 

- 

96.6 

69,0 

75.0 

286 

3067 

100 

55.4 

- 

• 

55.4 

61.1 

273 

3413 

9 

75.9 

- 

108 

75.9 

82.9 

323 

3799 

9 

89.0 

- 

116 

89.0 

97.7 

386 

4225 

9 

109 

- 

138 

109 

115 

575 

5000 

3 

207 

- 

- 

207 

187 

1160 

5000 

9 

126 

- 

155 

126 

130 

612 

5000 

100 

84.7 

- 

- 

84.7 

100 

480 

6000 

9 

182 

- 

- 

182 

196 

893 


TABLE II: Values of and (for 50% and 90% de- 
tection efficiencv), in units of 10“ 22 Hz -1 / 2 , for sine-Gaussian 
waveforms with the central frequency /o and quality factor 
Q. Three columns in the middle are the measured with 
the individual search algorithms for the H1H2L1V1 network. 
The next column is the of the logical OR of the cWB 
and II algorithms for the H1H2L1V1 network. The last two 
columns are the hfs S % and the of the logical OR of the 
algorithms and networks (H1H2L1V1 or H1H2L1 or H1H2). 
All h TB a values take into account statistical and systematic 
uncertainties as explained in Sec. VI. 
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/ T 

Hz] [trs| 

all networks, 
Lin. Circ. 

all networks, h 9 ^ 0 
Lin. Circ. 

1590 200 

34.7 

30.0 

131 

60.0 

2000 1.0 

49.5 

43.8 

155 

81.1 

2090 200 

43.3 

36.5 

155 

72.9 

2590 200 

58.6 

46.0 

229 

88.8 

3067 0.65 

88.2 

73.3 

369 

142 


TABLE III: Values of /&° % and h 9 ™ (f or 50 % and 90% de- 
tection efficiency using cWB), in units of 10~ 22 Hz~ 1/2 , for 
linearly and circularly polarized ringdowns characterized by 
parameters / and r. All h ras values take into account statis- 
tical and systematic uncertainties as explained in Sec. VI. 


/low 

A / 

r 

H 1 H 2 L 1 V 1 , 

all networks 

[Hz] 

[Hz] 

[ms] 

cWB 

a 

cWB or O 

h iO% 

'‘'res 

us? 

100 

100 

0.1 

7.6 

13.6 

7.6 

8.4 

19.6 

250 

100 

0.1 

9.1 

10.2 

8.8 

8.6 

18.7 

1000 

10 

0.1 

20.9 

28.6 

21.0 

21.8 

52.6 

1000 

1000 

0.01 

36.8 

38.2 

35.0 

36.3 

74.7 

1000 

1000 

0.1 

60.3 

81.7 

60.7 

63.5 

140 

2300 

100 

0.1 

40.4 

- 

40.4 

44.1 

94.4 

2000 

1000 

0.01 

60.7 

- 

60.7 

62.4 

128 

3500 

100 

0.1 

74.3 

- 

74.3 

84.8 

182 

3500 

1000 

0.01 

103 

- 

103 

109 

224 

5000 

100 

0.1 

101 

- 

101 

115 

255 

5 G 00 

1000 

0.01 

152 

- 

152 

144 

342 


TABLE IV: Values of /i“ % and (for 50% and 90% de- 
tection efficiency), in units of 10 -22 Hz~ 1/2 , for band-limited 
noise waveforms characterized by parameters /i ow , A /, and 
r. Two columns in the middle are the /i^ s % for the individ- 
ual search algorithms for the H1H2L1V1 network. The next 
column is the of the logical OR of the cWB and Q algo- 
rithms for the H1H2L1V1 network. The last two columns are 
the /ifs a % and the h 9 ^ A of the logical OR of the algorithms and 
networks (H1H2L1V1 or H1H2L1 or H1H2). All k rsa values 
take into account statistical and systematic uncertainties as 
explained in Sec. VI. 


higher h TSS by 11.1%. In the frequency band above 2 kHz, 
a new methodology, based on MonteCarlo simulations 
has been adopted to marginalize over calibration uncer- 
tainties: basically, we inject signals whose amplitude has 
been jittered according to the calibration uncertainties. 
The effect of miscalibration resulted in the increase of the 
combined h^° by 3 % to 14%, depending mainly on the 
central frequency of the injected signals. 


VII. SEARCH RESULTS 

In Section III we described the main steps in our search 
for gravitational-wave bursts. In the search all analysis 
cuts and thresholds are set in a blind way, using time- 


shifted (background) and simulation data. The blind cuts 
are set to yield a false-alarm rate of approximately 0.05 
events or less over the observation time of each search 
algorithm, network configuration and target frequency 
band. Here we describe the results. 


A. Candidate events 

After these cuts are fixed, the unshifted events are ex- 
amined and the various analysis cuts, DQFs, and vetoes 
are applied. Any surviving events are considered as can- 
didate gravitational-wave events and subject to further 
examination. The purpose of this additional step is to go 
beyond the binary decision of the initial cuts and evalu- 
ate additional information about the events which may 
reveal their origin. This ranges from “sanity checks” to 
deeper investigations on the background of the observa- 
tory, detector performances, environmental disturbances 
and candidate signal characteristics. 

Examining the unshifted data, we found one fore- 
ground event among all the different search algorithms 
and detector combinations that survives the blind se- 
lection cuts. It was produced by cWB during a time 
when all five detectors were operating simultaneously. As 
the possible first detection of a gravitational-wave signal, 
this event was examined in great detail according to our 
follow-up checklist. We found no evident problem with 
the instruments or data, and no environmental or instru- 
mental disturbance detected by the auxiliary channels. 
The event was detected at a frequency of 110 Hz, where 
all detectors are quite non-stationary, and where both 
the GEO 600 and Virgo detectors had poorer sensitivity 
(see Fig. 1). Therefore, while the event was found in the 
H1H2L1V1 analysis, we also re-analyzed the data using 
cWB and the H1H2L1 network. Figure 4 (top) shows the 
event above the blind selection cuts and the comparison 
with the measured H1H2L1 background of cWB in the 
frequency band below 200 Hz. 

No foreground event passes the blind selection cuts in 
the Cl H1H2L1 analysis (see Figure 4 (bottom)); more- 
over, there is no visible excess of foreground events with 
respect to the expected background. The cWB event is 
well within the tail of the Cl foreground and does not pass 
the final cut placed on correlated energy of the Hanford 
detectors. Furthermore, the event is outside of the fre- 
quency band (300-5000 Hz) processed by the EGC algo- 
rithm. Figure 5 (top) shows the corresponding EGC fore- 
ground and background distributions for the H1H2L1V1 
network. For comparison, Figure 5 (bottom) shows sim- 
ilar distributions from cWB, with no indication of any 
excess of events in the frequency band 1200-6000 Hz. 

To better estimate the significance of the surviving 
cWB event, we performed extensive background studies 
with cWB for the H1H2L1 network, accumulating a back- 
ground sample with effective observation time of approx- 
imately 500 years. These studies indicate an expected 
false alarm rate for similar events of once per 43 years 
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FIG. 4: Distribution of background (solid line) find fore- 
ground (solid dots) events from the search below 200 Hz in 
the H1H2L1 network, after application of category 2 data 
quality and vetoes: cWB (top), Q (bottom). The event- 
strength figures of merit on the horizontal axes are defined 
in the appendices on the search algorithms. The small error 
bars on the solid line are the 1 a statistical uncertainty on 
the estimated background, while the wider gray belt repre- 
sents the expected root-mean-square statistical fluctuations 
on the number of background events in the foreground sam- 
ple. The loudest foreground event on the top plot is the only 
event that survived the blind detection cuts of this search, 
shown as vertical dashed lines. This event was later revealed 
to have been a blind injection. 


for the cWB algorithm and the H1H2L1 network. The 
statistical significance of the event must take into account 
a “trials factor” arising from multiple analyses using dif- 
ferent search algorithms, networks and frequency bands. 
Neglecting a small correlation among the backgrounds, 
this factor can be estimated by considering the total effec- 
tive analyzed time of all the independent searches, which 
is 5.1 yr. The probability of observing one event at a 
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FIG. 5: Distribution of background (solid line) and fore- 
ground (solid dots) H1H2L1V1 events after category 2 data 
quality and vetoes: EGC events in the frequency band 300- 
5000 Hz (top), cWB events in the frequency band 1200- 
6000 Hz (bottom). The event-strength figures of merit on the 
horizontal axes are defined in the appendices on the search 
algorithms. The small error bars on the solid line are the 1 
a statistical uncertainty on the estimated background, while 
the wider gray belt represents the expected root-mean-square 
statistical fluctuations on the number of background events 
in the foreground sample. 


background rate of once per 43 years or less in any of our 
searches is then on the order of 10%. This probability was 
considered too high to exclude a possible accidental origin 
of this event, which was neither confirmed nor ruled out 
as a plausible GW signal. This event was later revealed 
to be a hardware injection with h rss = 1.0 x 10“ 2i Hz” 1 / 2 . 
It was the only burst injection within the “blind injection 
challenge.” Therefore it was removed from the analysis 
by the cleared injection data quality flag. We can report 
that cWB recovered the injection parameters and wave- 
forms faithfully, and the exercise of treating the event as 
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a real GW candidate was a valuable learning experience. 

Although no other outstanding foreground events were 
observed in the search, we have additionally examined 
events in the data set with relaxed selection cuts, namely, 
before applying category 3 DQFs and vetoes. In this set 
we find a total of three foreground events. One of these is 
produced by the EGC algorithm (0.16 expected from the 
background) and the other two are from the H-pipeline 
(1.4 expected). While an exceptionally strong event in 
the enlarged data set could, in principle, be judged to be 
a plausible GW signal, none of these additional events is 
particularly compelling. The EGC event occurred during 
a time of high seismic noise and while the H2 interferom- 
eter was re-acquiring lock (and thus could occasionally 
scatter light into the HI detector), both of which had 
been flagged as category 3 data quality conditions. The 
fl-pipeline events fail the category 3 vetoes due to having 
corresponding glitches in HI auxiliary channels. None of 
these three events passes the cWB selection cuts. For 
these reasons, we do not consider any of them to be a 
plausible gravitational-wave candidate. Also, since these 
events do not pass the predefined category 3 data qual- 
ity and vetoes, they do not affect the calculation of the 
upper limits presented below. 

B. Upper limits 

The S5y2/VSR1 search includes the analysis of eight 
network configurations with three different algorithms. 
We use the method presented in [27] to combine the 
results of this search, together with the S5yl search 
[12], to set frequentist upper limits on the rate of burst 
events. Of the S5y2 results, we include only the networks 
H1H2L1V1, H1H2L1 and H1H2, as the other networks 
have small observation times and their contribution to 
the upper limit would be marginal. Also, we decided a 
priori to use only the two algorithms which processed 
the data from the full S5y2 run, namely cWB and ft. 
(EGC only analyzed data during the ~5 months of the 
run when Virgo was operational.) We are left therefore 
with six analysis results to combine with the S5yl results 
to produce a single upper limit on the rate of GW bursts 
for each of the signal morphologies tested. 

As discussed in [27], the upper limit procedure com- 
bines the sets of surviving triggers according to which 
algorithm (s) and/or network detected any given trigger, 
and weights each trigger according to the detection effi- 
ciency of that algorithm and network combination. For 
the special case of no surviving events, the 90% confi- 
dence upper limit on the total event rate (assuming a 
Poisson distribution of astrophysical events) reduces to 
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FIG. 6: Selected exclusion diagrams showing the 90% con- 
fidence rate limit as a function of signal amplitude for Q=9 
sine-Gaussian (top) and Gaussian (bottom) waveforms for the 
results of the entire S5 and VSR1 runs (S5/VSR1) compared 
to the results reported previously (SI, S2, and S4). 


In the limit of strong signals in the frequency band be- 
low 2 kHz, the product Ct ot T is 224.0 day3 for S5yl and 
205.3 days for S5y2/VSR1. The combined rate limit 
for strong GW signals is thus 2.0 yr” 1 . For the search 
above 2 kHz, the rate limit for strong GW signals is 
2.2yr“ 1 . This slightly weaker limit is due to the fact 
that less data was analyzed in the S5yl high-frequency 
search than in the S5yllow-ffequency search (only 161.3 
days of H1H2L1 data [13]). Figure 6 shows the combined 
rate limit as a function of amplitude for selected Gaussian 
and sine-Gaussian waveforms. 

The results can also be interpreted as limits on the 
rate density (number per time per volume) of GWBs as- 
suming a standard-candle source. For example, given an 
isotropic distribution of sources with amplitude h rs& at 
a fiducial distance r 0 , and with rate density 7£, the rate 
of GWBs at the Earth with amplitudes in the interval 
[/i, h + dh ] is 


#90% — 


2.3 

ttotT ’ 


(7.1) 


dN = 


47r7£(/i r8S ro) 3 

h 4 


dh. 


(7.2) 


where 2.3 = log(l - 0.9), e tot is the detection efficiency 
of the union of all search algorithms and networks, and 
T is the total observation time of the analyzed data sets. 


(Here we have neglected the inclination angle i\ equiva- 
lently we can take h 2 to be averaged over cos l.) The ex- 
pected number of detections given the network efficiency 
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FIG. 7: Rate limit per unit volume at the 90% confidence 
level for a linearly polarized sine-Gaussian standard-candle 
with Egw = M 0 c 2 . 


e(/i) (for injections without any i dependence) and the 
observation time T is 


Adet 



= 4n7lT(h rss ra) 3 f 

Jo 


™dhh,- 4 e(h) 


(7.3) 


For linearly polarized signals distributed uniformly in 
cos t, the efficiency is the same with h rescaled by a fac- 
tor sin 2 l divi ded by that factor’s appropriately averaged 
value y/S/lE. Thus the above expression is multiplied 
by / (J 1 dcost(15/8) 3 / 2 sin 6 1 ^ 1.17. The lack of detec- 
tion candidates in the S5/VSR1 data set implies a 90% 
confidence upper limit on rate density 7Z of 


^ 90 % = 


47rT(/i rss 'r 0 ) 3 J^°dhh~ 4 e(h) ' 


Assuming that a standard-candle source emits waves 
with energy Egw — M©c 2 , where M© is the solar mass, 
the product h T3S r 0 is 



Figure 7 shows the rate density upper limits as a function 
of frequency. This result can be interpreted in the follow- 
ing way: given a source with a characteristic frequency / 
and energy Eq w — Me 2 , the corresponding rate limit 
is 7 ^oc%(/)(M©/M) 3 / 2 yr _1 Mpc“ 3 . For example, for 
sources emitting at 150 Hz with Eqw = 0.01M©c 2 , the 
rate limit is approximately 6 x l()" 4 yr~ 1 Mpc“ 3 . The 
bump at 361 Hz reflects the effect of the “violin modes” 
(resonant frequencies of the wires suspending the mir- 
rors) on the sensitivity of the detector. 


VIII. SUMMARY AND DISCUSSION 

In this paper we present results of new all-sky untrig- 
gered searches for gravitational wave bursts in data from 
the first Virgo science run (VSRl in 2007) and the second 
year of the fifth LIGO science run (S5y2 in 2006-2007). 
This data set represented the first long-term operation of 
a worldwide network of interferometers of similar perfor- 
mance at three different sites. Data quality and analysis 
algorithms have improved since similar searches of the 
previous LIGO run (S4 in 2004) [28] and even since the 
first year of S5 (S5yl in 2005-2006) [12, 13]. This is re- 
flected in an improved strain sensitivity with as low 
(good) as 5.6 x 10“ 22 Hz -1 / 2 for certain waveforms (see 
Table II), compared to best values of 1.3 x 10 -21 Hz -1 / 2 
and 6.0 x 10“ 22 Hz -1 / 2 for S4 and S5yl respectively. The 
new searches also cover an extended frequency band of 
50-6000 Hz. 

No plausible gravitational wave candidates have been 
identified in the S5y2/VSR1 searches. Combined with 
the S5yl results, which had comparable observation time, 
this yields an improved upper limit on the rate of bursts 
(with amplitudes a few times larger than A 3 g S % ) of 2.0 
events per year at 90% confidence for the 64-2048 Hz 
band, and 2.2 events per year for higher- frequency bursts 
up to 6 kHz. Thus the full S5/VSR1 upper limit is better 
than the S5yl upper limits of 3.75 per year (64-2000 Hz) 
and 5.4 per year (1-6 kHz), and is more than an order 
of magnitude better than the upper limit from S4 of 55 
events per year. 

We note that the IGEC network of resonant bar de- 
tectors set a slightly more stringent rate limit, 1.5 events 
per year at 95% confidence level [14]. However, those 
detectors were sensitive only around their resonant fre- 
quencies, near 900 Hz, and achieved that rate limit 
only for signal amplitudes (in h TS8 units) of a few times 
10 -19 Hz -ly/2 or greater, depending on the signal wave- 
form. (See Sec. X of [29] for a discussion of this com- 
parison.) Further IGEC observations during 6 months 
of 2005 [15] improved the rate limit to ~8.4 per year 
for bursts as weak as a few times 10~ 2O Hz -1//2 but did 
not change the more stringent rate limit for stronger 
bursts. The current LIGO-Virgo burst search is sensi- 
tive to bursts with A ras one to two orders of magnitude 
weaker than those which were accessible to the IGEC 
detectors. 

To characterize the astrophysical sensitivity achieved 
by the S5y2/VSR1 search, we calculate the amount of 
mass, converted into GW burst energy at a given dis- 
tance ro, that would be sufficient to be detected by the 
search with 50% efficiency (Mqw)- Inverting Eq. (7.5), 
we obtain a rough estimate assuming an average source 
inclination angle (i.e. A 2 SS is averaged over cost): 

Mgw = ^,o 2 /o 2 /4. (8.1) 

For example, consider a sine-Gaussian signal with /o = 
153 Hz and Q — 9, which (from Table II) has h™ % = 
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6.0 x 10 -22 Hz -1 / 2 for the four-detector network. As- 
suming a typical Galactic distance of lOkpc, that h T8S 
corresponds to Mow = 1.8 x 10 -8 M 0 . For a source in 
the Virgo galaxy cluster, approximately 16 Mpc away, the 
same would be produced by a mass conversion of 
roughly 0.046 A/ 0 . These figures are slightly better than 
for the S5yl search and a factor of ~5 better than the S4 
search. 

We also estimate in a similar manner a detection range 
for GW signals from core-collapse super novae and from 
neutron star collapse to a black hole. Such signals are 
expected to be produced at a much higher frequency 
(up to a few kHz) and also with a relatively small GW 
energy output (10“ 9 - 10“ 5 M 0 c ?). For a possible su- 
pernova scenario, we consider a numerical simulation of 
core collapse by Ott et al. [30]. For the model s25WW, 
which undergoes an acoustically driven explosion, as 
much as 8 x 10“ 5 M 0 may be converted to gravitational 
waves. The frequency content produced by this partic- 
ular model peaks around ~ 940 Hz and the duration is 
of order one second. Taking this to be similar to a high- 
er sine-Gaussian or a long-duration white noise burst, 
from our detection efficiency studies we estimate h^ Q 
of 17-22xl0 -22 Hz” 1 / 2 , i.e. that such a signal could be 
detected out to a distance of around 30 kpc. The ax- 
isymmetric neutron star collapse signals D1 and D4 of 
Baiotti et al. [5] have detection ranges (at 50% confi- 
dence) of only about 25 pc and 150 pc (see Fig, 3, due 
mainly to their lower energy (Mqw < 10” 8 Af 0 ) and 
also to emitting most of that energy at 2-6 kHz, where 
the detector noise is greater. 

The Advanced LIGO and Virgo detectors, currently 
under construction, will increase the detection range of 
the searches by an order of magnitude, therefore increas- 
ing by ~1000 the monitored volume of the universe. With 
that sensitivity, GW signals from binary mergers are 
expected to be detected regularly, and other plausible 
sources may also be explored. Searches for GW burst 
signals, capable of detecting unknown signal waveforms 
as well as known ones, will continue to play a central role 
as we increase our understanding of the universe using 
gravitational waves. 
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Appendix A: Data Quality Flags 

The removal of poor-quality LIGO data uses the data 
quality flag (DQF) strategy described in the first year 
analysis [12]. For the second year there are several new 
DQFs. New category 2 flags mark high currents in the 
end test-mass side coils, discontinuous output from a 
tidal compensation feed-forward system, periods when 
an optical table was insufficiently isolated from ground 
noise, and power fluctuations in lasers used to thermally 
control the radius of curvature of the input test masses. 
A flag for overflows of several of the main photodiode 
readout sensors that was used as a category 3 flag in the 
first year was promoted to category 2. New category 3 
flags mark noise transients from light scattered from HI 
into H2 and vice versa, large low-frequency seismic mo- 
tions, the optical table isolation problem noted above, 
periods when the roll mode of an interferometer optic 
was excited, problems with an optical level used for mir- 
ror alignment control, and one period when H2 was op- 
erating with degraded sensitivity. The total “dead time” 
(fraction of five time removed) during the second year of 
S5 due to category 1 DQFs was 2.4%, 1.4%, and <0.1% 
for HI, H2, and LI, respectively. Category 2 DQF dead 
time was 0.1%, 0.1%, and 0.6%, and category 3 DQF 
dead time was 4.5%, 5.5%, and 7.7%. Category 4 flags, 
used only as additional information for follow-ups of can- 
didate events (if any), typically flag one-time events iden- 
tified by Collaboration members on duty in the observa- 
tory control rooms, and thus are quite different between 
the first and second years. 

Virgo DQFs are defined by study of the general be- 
havior of the detector, daily reports from the control 
room, online calibration information, and the study of 
loud transient events generated online from the uncali- 
brated Virgo GW channel by the Qonline [31] program. 
Virgo DQFs include out-of-science mode, hardware injec- 
tion periods, and saturation of the current flowing in the 
coil drivers. Most of them concern a well identified de- 
tector or data acquisition problem, such as the laser fre- 
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quency stabilization process being off, photodiode satu- 
ration, calibration line dropouts, and loss of synchroniza- 
tion of the longitudinal and angular control. Some loud 
glitches and periods of higher glitch rate are found to be 
due to environmental conditions, such as increased seis- 
mic noise (wind, sea, and earthquakes), and 50 Hz power 
line ground glitches seen simultaneously in many mag- 
netic probes. In addition, a faulty piezo-electric driver 
used by the beam monitoring system generated glitches 
between 100 and 300 Hz, and a piezo controlling a mir- 
ror on a suspended bench whose cabling was not well 
matched caused glitches between 100 and 300 Hz and 
between 600 and 700 Hz. The total dead time in VSR1 
due to category 1 DQFs was 1.4%. Category 2 DQF dead 
time was 2.6%, and category 3 DQF dead time was 2.5% 
[32]. 


Appendix B: Event- by-event vetoes 

Event-by-event vetoes discard gravitational- wave 
channel noise events using information from the many 
environmental and interferometric auxiliary channels 
which measure non-GW degrees of freedom. Our 
procedure for identifying vetoes in S5y2 and VSR1 
follows that used in S5yl [12]. Both the GW channels 
and a large number of auxiliary channels are processed 
by the KleineWelle (KW) [33] algorithm, which looks 
for excess power transients. Events from the auxiliary 
channels which have a significant statistical correlation 
with the events in the corresponding GW channel are 
used to generate the veto time intervals. Candidate 
events identified by the search algorithms are rejected if 
they fall inside the veto time intervals. 

Veto conditions belong to one of two categories which 
follow the same notation used for data quality flags. Cat- 
egory 2 vetoes are a conservative set of vetoes targeting 
known electromagnetic and seismic disturbances at the 
LIGO and Virgo sites. These are identified by requiring 
a coincident observation of an environmental disturbance 
across several channels at a particular site. The resulting 
category 2 ‘data selection cuts are applied to all analyses 
described in this paper, and remove ~0.2% of analyz- 
able coincident live time. Category 3 vetoes make use 
of all available auxiliary channels shown not to respond 
to gravitational waves. An iterative tuning method is 
used to maximize the number of vetoed noise events in 
the gravitational- wave channel while removing a minimal 
amount of time from the analysis. The final veto list is 
applied to all analyses below 2048 Hz, removing ~2% of 
total analyzable coincident live time. 

An additional category 3 veto condition is applied to 
Virgo triggers, based on the ratio of the amplitude of an 
event as measured in the in-phase (P) and quadrature 
(Q) dark port demodulated signals. Since the Q channel 
should be insensitive to a GW signal, large Q/P ratio 
events are vetoed. This veto has been verified to be safe 
using hardware signal injections [34], with a loss of live 


time of only 0.036%. 

Appendix C: EGC burst search 

The Exponential Gaussian Correlator (EGC) pipeline 
is based on a matched filter using exponential Gaussian 
templates [35, 36], 

<&(*) = exp e 2 *^, (Cl) 

where /o is the central frequency and r 0 is the duration. 
Assuming that real GWBs are similar to sine-Gaussians, 
EGC cross-correlates the data with the templates, 

(c2> 

Here x(f) and $(/) are the Fourier transforms of the 
data and template, and S(f) is the two-sided noise power 
spectral density. N is a template normalization factor, 
defined as 



We tile the parameter space (/o, Qo = 2?rro/o) using the 
algorithm of [37]. The minimal match is 72%, while the 
average match between templates is 96%. The analysis 
covers frequencies from 300 Hz to 5 kHz, where LIGO 
and Virgo have comparable sensitivity. Qo varies from 2 
to 100, covering a large range of GW burst durations. 

The quantity p = y/2\C\ 2 is the signal-to-noise ra- 
tio (SNR), which we use to characterize the strength of 
triggers in the individual detectors. The analysis is per- 
formed on times when at least three of the four detec- 
tors were operating. Triggers are generated for each of 
the four detectors and kept if the SNR is above 5. In 
order to reduce the background, category 2 DQFs and 
vetoes are applied, followed by several other tests. First, 
triggers must be coincident in both time and frequency 
between a pair of detectors. The time coincidence win- 
dow is the light travel time between the interferometers 
plus a conservative 10 ms allowance for the EGC timing 
accuracy. The frequency coincidence window is selected 
to be 350 Hz. Second, events seen in coincidence in HI 
and H2 with a unexpected ratio in SNR are discarded 
(the SNR in HI should be approximately 2 times that in 
H2). Surviving coincident triggers are ranked according 
to the combined SNR, defined as 

Pcomb = yjp\ + P 2 j (C4) 

where pi and p 2 are the SNR in the two detectors. Third, 
a threshold is applied on p\ and p 2 to reduce the trigger 
rate in the noisier detector. This lowers the probability 
that a detector v/ith a large number of triggers will gen- 
erate manv coincidences with a few loud glitches in the 
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Network 

Obs. time 
[daysl 

# lags 

FAR 

Pcomb 

H1H2UV1 

66.6 

200 

< 400 Hz: 1 event in 10 years 

69.8 




> 400 Hz: 0.05 events 

21.0 

H1H2L1 

18.3 

1000 

< 400 Hz: 1 event in 10 years 

80.9 




> 400 Hz: 0.05 events 

10.0 

H1H2V1 

15.9 

1000 

< 400 Hz: 1 event in 10 years 

89.6 




> 400 Hz: 0.05 events 

15.4 

H1L1V1 

4.5 

2000 

< 400 Hz: 1 event in 10 years 

67.9 




> 400 Hz: 0.01 events 

24.2 


TABLE V: Thresholds and background tuning information 
for all the networks studied by the EGC pipeline. 


other detector. Finally, for each coincident trigger we 
compute the SNR disbalance measure 


A = 


Pcomb 

Pcomb + | Pi “ P 2 I 


(C5) 


This variable is useful in rejecting glitches in a pair of co- 
aligned detectors with similar sensitivity, and so is used 
primarily for pairs of triggers from the LIGO detectors. 

The background is estimated for each detector pair by 
time shifting the trigger lists. 200 time slides are done 
for H1H2L1V1, and more for the three-detector networks 
due tc their shorter observation times (see Table V). The 
thresholds applied to pi , P 2 and A are tuned for each de- 
tector pair to maximize the average detection efficiency 
for sine-Gaussian waveforms at a given false alarm rate. 
Once the pi, P 2 and A thresholds are applied, all trig- 
ger pars from the network are considered together and 
Pcomb is used as the final statistic to rank the triggers. A 
threshold is placed on p com b > chosen to give a low false 
alarm rate. More precisely, as we observe an excess of 
loud glitches with /o < 400 Hz, we use different thresh- 
olds depending on the frequency of the coincident trig- 
gers. Below 400 Hz, the false alarm rate is tuned to 1 
event per 10 years. Above, the threshold for each net- 
work is set to give a maximum of 0.05 events expected 
from background for that network. An exception is made 
for H1L1V1, where the maximum number is chosen to be 
0.01 events because of its shorter observation time. The 
final thresholds for each network are given in Table V. 


Appendix D: D-Pipeline burst search 

The H-Pipeline is essentially identical to QPipeline, 
which was used in previous LIGO S5 searches [12, 13]. 
QPipeline has since been integrated into a larger software 
suite, with a change in nomenclature but no significant 
change in methodology. Since this approach is discussed 
in detail in [12, 19], we provide only a summary here. 

The H-Pipeline, like EGC, functions as a matched- 
filter search on a single-interferometer basis. The data 
stream is whitened by linear predictive filtering [38] , then 
projected onto a template bank of complex exponen- 
tials. These templates are similar to those used by EGC, 


parametrized by central time tq, central frequency /o, 
and quality factor Q 0 > hut use bisquare windows rather 
than Gaussian windows. The template spacing is also 
different, selected for computational speed, rather than 
for strict mathematical optimization as in EGC. The tt 
template bank has a minimal match of 80%, and covers 
a frequency range from 48 Hz to 2048 Hz and a Q range 
from 2.35 to 100. 

The significance of a single-interferometer trigger is 
given by its normalized energy Z, defined as the ratio 
of the squared magnitude of X (the projection onto the 
best-matched template) for that trigger to the mean- 
squared magnitude of other templates with the same /o 
and Qo- Eor Gaussian white noise, Z is exponentially 
distributed and related to the matched filter SNR p by 

Z = \X\ 2 /(\Xf)=p>/2. (Dl) 

Z is used to rank LI and VI triggers. 

For Hi and H2, fi-Pipeline takes advantage of their 
co-located nature to form two linear combinations of the 
data streams. The first of these, the coherent stream 
H+, is the sum of the strains in the two interferometers 
weighted by their noise power spectral densities. We de- 
fine the coherent energy following (Dl). We also 
define the correlated energy Zg”* r , which is obtained by 
removing the contribution to Z^ + from HI and H2 indi- 
vidually and leaving only the cross-correlation term[42]. 
The H1H2 cuts are based on Z^ r , because it is less sus- 
ceptible than Zh + to instrumental glitches, so providing 
better separation between signal and noise. The second 
stream, the null stream H_, is the difference between 
the strains in HI and H2. The normalized energy 
should be small for a gravitational wave, but generally 
much larger for an instrumental glitch. We therefore veto 
coherent stream triggers which are coincident in time and 
frequency with null stream triggers. 

We require triggers to be coincident in at least two 
detectors. The interferometer combinations analyzed are 
shown in Table VI. (Note that because of the coherent 
analysis of HI and H2, both must be operating for data 
from either to be analyzed.) Triggers are required to be 
coincident in both time and frequency as follows: 

|Ti -T 2 1 < T c + 1 max(<7i, <r 2 ) (D2) 

IfV-Fal < (D3) 

Here T and F are the central time and frequency of the 
triggers, a and b are their duration and bandwidth, and 
T c is the light travel time between the interferometers. 

The background for each detector pair is determined 
by time-shifting the triggers from one detector. We use 
1000 shifts for each pair, except H1-H2. Only 10 shifts 
between HI and H2 are used because the coherent anal- 
ysis requires each shift to be processed independently, 
substantially increasing the computational cost. Also, 
time shifts between HI and H2 are less reliable because 
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they miss correlated background noise from the common 
environment. For all pairs, triggers below and above 200 
Hz are treated separately because of the different char- 
acteristics of the glitch populations at these frequencies. 

Normalized energy thresholds are set separately for 
each detector combination and frequency range such that 
there is less than a 5% probability of a false alarm after 
category 3 DQFs and vetoes. Table VI shows the thresh- 
olds fcnd surviving events in timeslides for each combina- 
tion. Fig. 8 shows background and injection triggers and 
the energy thresholds for one interferometer pair. 


Detector combination 
and frequency band 

threshold 

events in 1000 
timeslides 

H1E2L1 < 200 Hz 

Zg 5 " > 37, Zli > 13 

14 

H1H2L1 > 200 Hz 

ZTJ > 13, Zli > 13 

16 

H1H2V1 < 200 Hz 

Zg° + rr > 22, Zvx > 13 

9 

H1H2V1 > 200 Hz 

Zg°; r > 14, Zvx > 13 

0 

L1V1 < 200 Hz 

Zli > 32 and Zvi > 
(4.7 x 10 _13 Zli)“ 0,3 

4 

L1V1 > 200 Hz 

Zli > 30 and Zvi > 
(6.9 x 10" 12 Zz,ir o>27 

5 

H1H2 < 200 Hz 

>80 

0 (10 slides) 

H1H2 > 200 Hz 

z§z r > 30 

0 (10 slides) 

total events 


48 


TABLE VI: Thresholds on normalized energy for the various 
detector combinations. 



Injections 
- Background 


FIG. 8: Distribution of background and injection triggers be- 
low 200 Hz after category 3 DQFs and vetoes for Ll-Vl pair. 
The dashed lines show the final normalized energy thresholds. 


Appendix E: Coherent WaveBurst Search 

Coherent WaveBurst (cWB) is a coherent algorithm 
for detecting gravitational-wave bursts. It constructs a 


least-squares fit of the two GW polarizations to the data 
from the different detectors using the constraint likeli- 
hood method [25]. The cWB algorithm was first used in 
search for gravitational wave bursts in the LIGO-GEO 
network [39]. More recently it has been used in the LIGO 
S5 first-year low-frequency search [12], and detailed de- 
scriptions of the algorithm can be found there and in 
[24, 25]. 

The cWB analysis in this search covers frequencies 
from 64 Hz to 6.0 kHz, with the data processing split 
into two bands. The low-frequency (LF) band (64 Hz to 
2.0 kHz) contains the most sensitive (but also the most 
non-stationary) data. The high-frequency (HF) band 
(1.28 kHz to 6.0 kHz) is dominated by the shot noise of 
the detectors and is much less polluted by environmen- 
tal and instrumental transients. Splitting the analysis 
into two bands is convenient for addressing the differ- 
ent noise characteristics in these bands. It also eases the 
computational cost. The overlap of the bands is used to 
cross-check the results and to preserve the sensitivity to 
wide-band signals near the boundary between the bands. 

The cWB analysis is performed in several steps. First, 
the data are decomposed into Meyer wavelets. Time- 
frequency resolutions of (8 x 1/16. 16 x 1/32, 32 x 1/64, 
64 x 1/128, 128 x 1/256, 256 x 1/512 [Hz x s]) are used 
for the low-frequency search and (12.5 x 1/25, 25 x 1/50, 
50 x 1/100, 100 x 1/200, 200 x 1/400, 400 x 1/800 [Hz 
x s]) for the high-frequency search. The data are pro- 
cessed with a linear predictor error filter to remove power 
lines, violin modes and other predictable data compo- 
nents. Triggers are identified as sets of wavelet pixels 
among the detectors containing excess power at time de- 
lays consistent with a gravitational wave from a phys- 
ical sky position. For each trigger, trial incoming sky 
locations are sampled with 1° resolution, and various co- 
herent statistics are computed. These include the maxi- 
mum likelihood ratio L m (a measure of the sum-squared 
matched-filter SNR detected in the network), the net- 
work correlated amplitude 77, the network correlation co- 
efficient (cc), the energy disbalance statistics A, Ahh and 
the penalty factor Pf. (Each of these statistics is de- 
scribed in detail in [12].) The trial sky position giving 
the largest cc • is selected as the best-guess incident 
direction, and the coherent statistics for this position are 
recorded. Finally, several post-production selection cuts 
are applied to the triggers to reduce the background. 

Two groups of selection cuts are used in cWB. First, 
cuts on cc, A, Ahh and Pf are used to distinguish noise 
outliers from genuine GW signals. The most powerful 
consistency cut is based on the network correlation co- 
efficient cc. For example, Figure 9 shows a scatter plot 
of background triggers as a function of 77 and cc. Strong 
outliers (with large values of 77) are characterized by low 
values of cc and are well separated from simulated sig- 
nals. Additional selection cuts are based on the energy 
disbalance statistics Anet> Ahh, and Pf. They are used 
to reject specific types of background events, such as Hl- 
H2 correlated transients. 
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The second, final cut is on the network correlated 
amplitude 77 , which characterizes the significance of the 
triggers. Due to different characteristics of the back- 
ground noise during the run and in the different fre- 
quency bands, the threshold on 77 is selected separately 
for each network configuration and frequency band to 
give false alarm probabilities of a few percent. In partic- 
ular, in the low-frequency search separate 77 thresholds 
are used for triggers below and above 200 Hz. Table VII 
shows the thresholds used in the analysis. 


FIG. 9: Distribution of background triggers (black dots) after 
category 2 DQFs and vetoes for the L1H1H2 network for the 
high-frequency search, with Q = 9 sine-Gaussians injections 
(gray dots). The dashed lines show the thresholds on r\ and 
cc chosen for this network. 


The background is estimated separately on each seg- 
ment of data processed. The cWB algorithm forms cir- 
cular data buffers and shifts one detector with respect to 
the others, repeating the analysis hundreds of times on 
the time-shifted data. Table VIII shows the number of 
lags and accumulated background observation time for 
the various cWB searches and network configurations. 
The background data sets are used for tuning of the cWB 
selection cuts and also for estimation of the significance 
of the foreground events. For example, to estimate the 
significance of the blind injection identified by cWB, we 
generated a background sample with observation time 
equivalent to approximately 1000 HIH2L1 S5y2/VSR1 
data sets. 
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cWB LF 

cc 771 772 Ahh Anet Pj 

cWB HF 

CC 77 Pj 

H1H2L1V1 

0.5 4.5 4.1 

0.3 

0.35 

0.6 

0.6 4.3 0.6 

H1H2L1 

0.6 6.0 4.2 

0.3 

0.35 

0.6 

0.6 4.1 0.6 

II1H2V1 

- - - 

- 

- 

- 

0.7 4.6 0.6 

H1L1V1 

0.5 5.0 5.0 

- 

0.7 

0.6 

0.6 6.0 0.6 

H1H2 

0.6 6.0 4.1 

0.3 

0.35 

0.6 

0.5 5.0 0.6 

HILl 

0.6 9.0 5.5 

- 

0.35 

0.6 

- - - 

H2L1 

0.6 6.5 5.5 

- 

0.35 

0.6 

- - - 


TABLE VII: Thresholds for each network for the cWB low- 
and high-frequency searches. Different thresholds on q are 
used for triggers below 200 Hz (771) and above 200 Hz (772) 
due to the different characteristics of the LIGO background 
noise in these frequency ranges. No energy disbalance cuts 
(Anet> Ahh) are applied in the high-frequency search. In 
addition, a penalty factor cut of Pj > 0.6 is applied to all 
network configurations and searches. 
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cWB LF Background 

cWB HF Background 


# lags 

Obs time 

FAR 

#lags 

Obs time 

FAR 



[year] 

[year -1 ] 


[year] 

[year -1 ] 

H1H2L1V1 



0.3 

96 

17.6 

0.17 

H1H2L1 

1000 

499.9 

0.1 

96 

31.7 

0.09 

H1H2V1 

- 

- 

- 

96 

3.7 

0.27 

H1L1V1 

200 

1.8 

3.3 

288 

3.0 

0.33 


200 

28.2 

0.04 

192 

17.2 

0.06 


200 

4.6 

2.0 

- 

- 

- 


200 

1.6 

0.6 

- 

- 

- 


TABLE VIII: Background observation time and false alarm 
rates for each network for the cWB low- and high-frequency 
searches. 
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